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We propose a theoretical framework for analyzing two-body nonleptonic D meson decays, based 

on the factorization of short-distance (long-distance) dynamics into Wilson coefficients (hadronic 

matrix elements of four-fermion operators). The parametrization of hadronic matrix elements in 

terms of several nonperturbative quantities is demonstrated for the D — >■ PP decays, P denoting 

a pseudoscalar meson. We consider the evolution of Wilson coefficients with energy release in 

fNj individual decay modes, and the Glauber strong phase associated with the pion in nonfactorizable 

,-H annihilation amplitudes, that is attributed to the unique role of the pion as a Nambu-Goldstone 

(*^ boson and a quark-anti-quark bound state simultaneously. The above inputs improve the global 

fN| fit to the branching ratios involving the rj' meson, and resolves the long-standing puzzle from the 

^m. D" — > TT^Tv~ and D" — > K'^K~ branching ratios, respectively. Combining short-distance dynamics 

ir^ associated with penguin operators and the hadronic parameters determined from the global fit 

r^ to branching ratios, we predict direct CP asymmetries, to which the quark loops and the scalar 

penguin annihilation give dominant contributions. In particular, we predict AAcp = Acp{K^K~) — 

■^ Acp(7r+7r") = -1.00 X 10"^, lower than the LHCb and CDF data. 

CN 

^ I. INTRODUCTION 

a, 

O , Two-body nonleptonic D meson decays have attracted great attention recently. The very different data for the 

(]J D^ — > vr+TT^ and D^ — > K^K^ branching ratios have been a long-standing puzzle: the former (latter) is lower 

T^ (higher) than theoretical predictions from analyses based on the topology parametrization [TJ [5] . The deviation in 

these two D -^ PP modes, P representing a pseudoscalar meson, stands even after taking into account flavor SU{3) 

rf-s symmetry breaking effects in emission amplitudes [I]. Another puzzle appears in the difference between the direct 

J> CP asymmetries of the D^ ^- K^K^ and D^ — !• tt+tt" decays, 

o 

CN A^CP = Acp{K+K-) - Acp(7r+7r") = [-0.82 ± 0.21(stat) ± 0.11(syst)]%, (1) 

^H 

C^ measured by LHCb collaboration 3 , which is the first evidence of CP violation in charmed meson decays. The above 
(y~j result has been confirmed by the CDF collaboration |^^, 

o 

(N AAcP = [-0.62 ± 0.21(stat) ± 0.10(syst)]%. (2) 



X 



The quantity A^cp is naively expected to be much smaller in the Standard Model (SM), since the responsible penguin 
contributions are suppressed by both the CKM matrix elements and the Wilson coefficients [6l [^ , 



H A iKbKfcl as 4 .o^ 

The above two puzzles have triggered intensive investigations in the literature employing different approaches. 
The mechanism responsible for the very different D^ -^ tt+tt" and D" — >■ K'^K~ branching ratios is still unclear. 
It has been speculated that the long-distance resonant contribution through the nearby resonance /o(1710) in the 
W-exchange topology could explain why a 13° meson decays more abundantly to the K^K~ than tt+tt" final state 
[T]: it is attributed to the dominance of the scalar glueball content in /o(1710) and the chiral suppression on the 
scalar glueball decay into two pseudoscalar mesons. This mechanism was systematically formulated in the single pole- 
dominance model with additional arbitrary "Wilson coefficients" , and then included into the topology parametrization 
for global fits [8j. Unfortunately, the predicted D° — >■ tt+tt" (D° — >• K^K~) branching ratio remains higher (lower) 



1 The CDF collaboration observed the smaller direct CP asymmetries Acp(-D° -^ tt+tt-) = (+0.22 ± 0.24 ± 0.11)% and Acp(D° 
K+K-) = (-0.24 ± 0.22 ± 0.09)% in 0. 



than the data. That is, a systematic and global understanding of the D — )■ PP decays has not existed yet. Because 
of the high precision of the two data, even a small deviation deteriorates the global fit. If the branching ratios are 
not well explained, some important decay mechanism may be still missing, and any predictions for the direct CP 
asymmetries are not convincing. 

For the second puzzle, a reliable evaluation of the penguin contributions to two-body nonleptonic D meson decays 
is not available so far. In [9l [10] the tree amplitudes were determined by fitting the topology parametrization to 
the data of the branching ratios, while the penguin amplitudes were calculated in the QCD-improved factorization 
(QCDF) |TT1[T2]. That is, the tree and penguin contributions were treated in the different theoretical frameworks. It 
has been noticed that the penguin amplitudes derived from QCDF lead to a tiny AAqp of order 10~^ [TU]. Allowing 
the penguin amplitudes to be of the same order as the tree ones discretionally, AAcp reaches —0.13% ~ O(10~^) 
[lOj . In another work |13j also based on the topology parametrization, the penguin contribution via an internal b 
quark was identified as the major source of CP violation, since it cannot be related to the tree amplitudes. This 
penguin contribution, including its strong phase, was constrained by the LHCb data, and then adopted to predict 
direct CP asymmetries of other decay modes. In this way it is difficult to tell whether the large A^cp ~ O(10~^) in 
Eqs. (fTl) and (21) arise from new physics [TH - ED] . The similar comment applies to the analysis in flV. Viewing that 
the postulated AAqp ranges from 10~^ to 10~^ in the SM [22,, it is crucial to give precise predictions for the direct 
CP asymmetries in the D -^ PP decays. 

In this paper we aim at proposing a theoretical framework, in which various contributions, such as emission and 
annihilation ones, and tree and penguin ones are all handled consistently. The only assumption involved is the 
factorization of short-distance dynamics and long-distance dynamics in two-body nonleptonic D meson decays into 
Wilson coefficients and hadronic matrix elements of four-fermion operators, respectively. For the hadronic matrix 
elements, including the emission, M^-annihilation, and M^-exchange amplitudes, the treatment is similar to the topology 
parametrization in the factorization hypothesis. The quark-loop and magnetic penguin contributions are absorbed 
into the Wilson coefficients for the penguin operators. As to the scale of the Wilson coefficients, we set it to the energy 
release in individual decay modes, which depends on masses of final states. An important ingredient is the Glauber 
strong phase factor [53] associated with a pion, whose strength has been constrained by the data of the B — >■ ttK 
direct CP asymmetries [24 . Briefly speaking, we have improved the topology parametrization by taking into account 
mode-dependent QCD dynamics, for instance, flavor SU{5) symmetry breaking effects, in two-body nonleptonic D 
meson decays. 

Due to the small charm quark mass nic ~ 1.3 GeV, a perturbative theory may not be valid for D meson decays. 
However, with the abundant data, all the above hadronic parameters can be determined. It will be shown that the 
Glauber phase associated with a pion is crucial for understanding the D^ — > tt+tt" branching ratio very different 
from the Z?° — >■ K'^K" one. This additional phase modifies the relative angle and the interference between the 
emission and annihilation amplitudes involving pions. The predicted D^ — >■ tt^tt^ branching ratio is then reduced, 
while the predicted £>+ — > tt+t] branching ratio is enhanced, such that the overall agreement with the data is greatly 
improved. Our scale choice for the Wilson coefficients impacts the D — > PP decays containing the t]' meson in the 
final states: the rj' meson mass about 1 GeV is not negligible compared to the D meson mass, so the energy release in 
these modes should be lower. Once the hadronic parameters are fixed by the measured branching ratios, the penguin 
amplitudes, expressed as the combination of the hadronic parameters and the corresponding Wilson coefficients, are 
also fixed in our framework. That is, we can predict the direct CP asymmetries of the D — )■ PP decays in the SM 
without ambiguity. In particular, we obtain AAqp — —1.00 x 10~^, which discriminates the opposite postulations 
on large (small) direct CP asymmetries in singly Cabibbo-suppressed D meson decays [5S] ([IS]). It will be observed 
that the quark loops and the scalar penguin annihilation generate large relative phases between the tree and penguin 
amplitudes, and are the dominant sources of the direct CP asymmetries. 

In Sec. II we present our parametrization of the tree contributions to the D — >■ PP branching ratios based on 
the factorization hypothesis with QCD dynamics being implemented. In Sec. Ill the penguin contributions from the 
operators O3-6, from O1.2 through the quark loops, and from the magnetic penguin Ogg are formulated. The direct 
CP asymmetries in the D — >■ PP decays are then predicted. Section IV is the conclusion. We show the explicit 
topology parametrization for the D" -^ tt+tt" and D*^ — > K^K^ decays in Appendix |AJ Appendix \B\ collects the 
evolution formulas for the Wilson coefficients in the scale fi. The 77-77' mixing formalism and the chiral factors of their 
fiavor states are given in Appendix [C] The form factors involved in the factorizable scalar penguin annihilation are 
specified in Appendix JP] 
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FIG. 1: Four topological diagrams contributing to D — >■ PP decays: (a) the color-favored tree amplitude T, (b) the color- 
suppressed tree amplitude C, (c) the W-exchange amplitude E, and (d) the W^-annihilation amplitude A. The thick line 
represents the charm quark. 



II. BRANCHING RATIOS 



In this section we formulate the tree contributions, which dominate the branching ratios of charm decays. The 
relevant weak effective Hamiltonian is written as 



n 



eff 



Gp 

75 



VcKM [Ci{ii)Oi{ii) + C2{ii)02{y)] + h.c, 



(4) 



where Gp is the Fermi coupling constant, Vckm denotes the corresponding Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements, Ci^2 are the Wilson coefficients, and the current-current operators are 



01 = iUaq2l3)v-Aiqil3Ca)v~A, 

02 = iUaq2a)v-A{qif3C/3)v~A, 



(5) 



with a, /? being color indices, and gi.2 being the dor s quark. There are four types of topologies dominantly contributing 
to the branching ratios of _D — > P1P2 decays [27], where P2 represents the meson emitted from the weak vertex: 
the color-favored tree amplitude T, the color-suppressed amplitude C, the PF-exchange amplitude E, and the W- 
annihilation amplitude A, shown in Fig. [T] Penguin contributions are neglected for branching ratios due to the 
suppression of the CKM matrix elements V*fjVub- Besides, there exist flavor-singlet amplitudes, in which a quark- 
antiquark pair produced from vacuum generates color- and flavor-singlet states like 77'''. However, they are also 
neglected in this work, since the current data do not suggest significant flavor-singlet contributions [J [2] . 

It is known that the heavy quark expansion does not work well for charm decays because of the small rric, so that 
nonfactorizable contributions, such as vertex corrections, final-state interaction (FSI), resonance effects, etc, cannot 
be ignored, especially in the amplitude C. We parameterize these contributions into a factor Xnf ■ For simplicity, 
we drop nonfactorizable contributions in the amplitude T, which is dominated by factorizable ones. In general, there 
exists a relative strong phase (j) between T and C, which may arise from inelastic FSI or other long-distance dynamics 
in D meson decays ^ [5] . Therefore, the emission amplitudes T and G are parameterized as 



{PlP2\'Heff\D)T. 



Gf 

V2 



VcKMaiMf^) fp^imj) - mpjF^^^{mp^ 



(6) 



with the Pi (P2) meson mass mp-^ (mp^), the P2 meson decay constant /p^, and the D -^ Pi transition form factor 



Fq ^ . The scale-dependent Wilson coefficients are given by 



1 



Xnfe'^ 



(7) 



for T and C, respectively, with Nc — 3 being the number of colors. 

The flavor SU{3) symmetry breaking effects have been known to be significant in the D — > PP decays, especially 
in the singly Cabibbo-suppressed modes [TJ |H1 |2H] ■ In additional to the different decay constants fp^ and form factors 
Fq \ the mass ratios mj^„(,) /niD should be distinguished too, with niK, tti^{,), and mD being the masses of the 
kaon, the r^^'^ meson, and the D meson, respectively. As suggested in the perturbative QCD (PQCD) approach |29j . 
the scale fj, is set to the energy release in individual decay processes, which depends on final-state masses. We propose 
the choice 



/^=^Amz3(l-r|), (8) 

where r2 = rri^p l"rr?T:) is the mass ratio of the P^ meson over the D meson, and A, representing the momentum of soft 
degrees of freedom in the D meson, is a free parameter. The evolution of the Wilson coefficients in the scale /i is 
shown in AppendixlBJ The above choice further takes into account the 5C/(3) breaking effect in the Wilson coefficients 
for the emission amplitudes. It is then legitimate to regard Xn/7 i> and A as being universal to all modes, which will 
be determined by experimental data. 

The D^ -^ K^K branching ratio through the pure T4^-exchange channel vanishes in the SU{3) limit due to the 
cancelation of the CKM matrix elements [U [30] . Hence, the large observed branching ratio of this mode manifests 
the breaking of the ffavor SU(3) symmetry in the VF-exchange contribution, which becomes almost of the same order 
as the emission one. Both the V7-exchange topology for neutral D meson decays and the W-annihilation topology for 
charged D meson decays will be included in our framework. Since the factorizable contributions to these amplitudes 
are down by the helicity suppression, we consider only the nonfactorizable contributions. Then the former (latter) is 
proportional to the Wilson coefficient C2 (C'l). Based on the above reasoning, we parameterize the amplitudes E and 
A as 

{P^P2\Ueff\D)E,A = ^VcKMhlf{ti)fDm], ( '-^^ ] , (9) 

with the factors 



= ^VcKMh'^f{n)fDml 


(fpjp. 

\ fi 


C(a^) = ^2(/i)xf.e<s 




fe^,(M) = Ci(A.)x^,e<=. 





(10) 

The dimensionless parameters Xqs ^^^^ ^f'ij describe the strengths and the strong phases of the involved matrix 
elements, where the superscripts g, s differentiate the light quarks and the strange quarks strongly produced in pairs. 
The &'s are defined for the tttt final states, to which those for other final states are related via the ratio of the decay 
constants fp^fp^/ f"^- That is, the flavor SU{2>) breaking effects from the strongly produced quark pairs and from the 
decay constants have been taken into account. The parameters x^s and (t>^'g are then assumed to be universal, and 
will be determined from data. Similarly, we consider the S'C/(3) breaking effect in the energy release, which defines 
the scale of the Wilson coefficients Ci_2(m) 

/i = ^KmD{l-rl){l-rl), (11) 

with the mass ratios ri_2 = ^np^^p^/mp,. 

As pointed out in [21], a special type of soft gluons, called the Glauber gluons, may introduce an additional strong 
phase to a nonfactorizable amplitude in two-body heavy-meson decays. The Glauber phase associated with a pion 
is more significant, because of its simultaneous role as a Nambu-Goldstone boson and a qq bound state: the valence 
quark and anti-quark of the pion are separated by a short distance in order to reduce the confinement potential 
energy, while the multi-parton states of the pion spread over a huge space-time in order to meet the role of a massless 
Nambu-Goldstone boson [31]. The relatively larger soft effect from the multi-parton states in the pion results in the 
Glauber phase. Hence, we assign a phase factor exp(iS'Tr) for each pion involved in the nonfactorizable amplitudes 



TABLE I: Branching ratios in units of percentages for Cabibbo-favored D — >■ PP decays. Our results are compared to those from 
the analysis including FSI effects [26], the topological-diagram approach [l], the pole-dominant model [8], and the experimental 
data 1321. 



Modes 


Br(FSI) 


Br(diagram) 


Br(pole) 


Br(exp) 


Br(this work) 


23° ^ ^«Z" 


1.35 


2.36±0.08 


2.4 ±0.7 


2.38±0.09 


2.41 


75« ^ ■,+K- 


4.03 


3.91±0.17 


3.9±1.0 


3.891±0.077 


3.70 


D° ^ k\ 


0.80 


0.98±0.05 


0.8 ±0.2 


0.96±0.06 


1.00 


Z3« ^ k\' 


1.51 


1.91±0.09 


1.9 ±0.3 


1.90±0.11 


1.73 


D+ -> 7r+F° 


2.51 


3.08±0.36 


3.1 ±2.0 


3.074±0.096 


3.22 


D+ ^ 7^+Z° 


4.79 


2.97±0.32 


3.0 ±0.9 


2.98±0.08 


3.00 


■Ds -^ "^^"n 


1.33 


1.82±0.32 


1.9±0.5 


1.84±0.15 


1.65 


D+ -^ TT+ri' 


5.89 


3.82±0.36 


4.6 ±0.6 


3.95±0.34 


3.44 


D+ ^ tt+tt" 










<0.06 






E and A. The phase factor associated with the tttt final states is then given by exp(2iSTr) as explicitly shown in 
Appendix lAj The reason why we do not introduce the Glauber phase into the nonfactorizable emission diagrams, but 
only into the VF-exchange and Vl^- annihilation diagrams is that the emission topology mainly receives factorizable 
contributions in most of the decay modes. Since the Glauber phase appears only in nonfactorizable amplitudes, we 
will not include it into the emission topology in the present work. 

In summary, our parametrization for the emission amplitudes in Eq. mf and the VF-exchange and VF-annihilation 
ones in Eq. ^ contains 12 free parameters. The global fit to all the data of 28 13 — >■ PP branching ratios leads to 
the outcomes 



A = 0.56 GeV, 
S^ = -0.50, 



Xnf = -0.59 
b = -0.62, 6 



= 4. 



= 0.11, 



xf = 0.18, 



= 4.23, 



X^ = 0.12, xf-0.17, 
= 4.06, (jjf = 3.48, 



(12) 



with the fitted x^ = 6.9 per degree of freedom, which marks a great improvement compared to x^ = 87 in [Ij. The 
parameter A = 0.56 GeV indeed corresponds to the soft momentum involved in the D meson, implying the typical 
energy release /i- '--^ 1 GeV. This scale, which is not very low, indicates that the PQCD factorization formulas [34] 
might provide a useful guideline for formulating various contributions to the D — ^ PP decays in our framework. Since 
the fermion flows of the amplitudes E and A can be converted into each other through the Fierz transformation, the 
relations Xq ~ xf: xf ~ xfj ^-^d the similar relations for the strong phases are expected. The satisfactory fit also 
confirms our postulation that the helicity suppression works well, and the nonfactorizable contributions dominate 
the W^-exchange and W^-annihilation amplitudes. A large SU{3) breaking effect has been revealed by the inequality 
Xq ¥" Xs ■ The values in Eq. (12) support i? > A in [TJ [5] according to Eq. (10 1, consistent with the fact that 
the D° meson lifetime is shorter than the _D+ meson lifetime. The Glauber phase Stt = —0.50 is in agreement with 
that extracted in [21], which resolves the puzzle from the dramatically different data for the direct CP asymmetries 
Acp{B° -^ tt'^K^) and Acp{B^ -^ tt^K^). 

Our results for the D — > PP branching ratios are presented in Tables [ij [llj and |III| for the Cabibbo-favored, singly 
Cabibbo-suppressed, and doubly Cabibbo-suppressed decays, respectively. The experimental data, as well as the 
predictions from other approaches, such as the topological-diagram approach [Tj, the factorization approach combined 
with the pole-dominant model for the M^-exchange and M^-annihilation contributions [8], and an analysis with FSI 
effects from nearby resonances [25j, are also listed for comparison. It is obvious that the predictions for the singly 

Cabibbo-suppressed decays, including the D^ — )■ tt^tt~ ^ K^K~ ^ and K'^K modes, are well consistent with the data: 
the known large S'C/(3) breaking effects in these decays have been properly described in our parametrization. The 



major sources to x^ of our fit come from the D" -^ tt'^tt" (Ax^ = 23) and -D+ 



(Ax-^ = 17) decays. These 



modes are exceptional in the sense that they involve the colo-suppressed emission amplitudes C, which are dominated 
by the nonfactorizable contributions. To reduce x^i more SU{?>) breaking effects, such as the Glauber phase in C 
[24], need to be introduced. 
The Wilson coefficients for the D'^ -^ tt+tt^ and K^K~ decays. 



ai(7r7r) = 1.09, 
ai{KK) = 1.10, 



^il47.8° 



a2(7r7r) =0.81e* 
a2(ifi^) = 0.83e*"«-2°. 



(13) 



TABLE II: Same as Table kl except for singly Cabibbo-suppressed D — >■ PP decays in units of 10 



Modes 


Br(FSI) 


Br (diagram) 


Br(pole) 


Br(cxp) 


Br(this work) 


D" ^ 7r+^- 


1.59 


2.24±0.10 


2.2 ±0.5 


1.45±0.05 


1.43 


D° ^ K+K~ 


4.56 


1.92±0.08 


3.0 ±0.8 


4.07±0.10 


4.19 


DO ^ K°K° 


0.93 





0.3 ±0.1 


0.320±0.038 


0.36 


D° ^ tv^tt" 


1.16 


1.35±0.05 


0.8 ±0.2 


0.81±0.05 


0.57 


23° ^ 7T% 


0.58 


0.75±0.02 


1.1 ±0.3 


0.68±0.07 


0.94 


D° ^ TfOj?' 


1.7 


0.74±0.02 


0.6 ±0.2 


0.91±0.13 


0.65 


D° ~^r)rj 


1.0 


1.44±0.08 


1.3 ±0.4 


1.67±0.18 


1.48 


D° -^ 7]T]' 


2.2 


1.19±0.07 


1.1±0.1 


1.05±0.26 


1.54 


D+ ^ n+Tv" 


1.7 


0.88±0.10 


1.0 ±0.5 


1.18±0.07 


0.89 


D+ -> K+k" 


8.6 


5.46±0.53 


8.4 ±1.6 


6.12±0.22 


5.95 


D+ -> 7r+J7 


3.6 


1.48±0.26 


1.6 ±1.0 


3.54±0.21 


3.39 


D+ ^ -K+ri' 


7.9 


3.70±0.37 


5.5 ±0.8 


4.68±0.29 


4.58 


D+ -> 7T°K + 


1.6 


0.86±0.09 


0.5 ±0.2 


0.62±0.23 


0.67 


D+ -> 7v+K° 


4.3 


2.73±0.26 


2.8 ±0.6 


2.52±0.27 


2.21 


D+ -^ K+r, 


2.7 


0.78±0.09 


0.8 ±0.5 


1.76±0.36 


1.00 


D+ -^ K+rj' 


5.2 


1.07±0.17 


1.4 ±0.4 


1.8±0.5 


1.92 



TABLE III: Same as Table ^1 except for doubly Cabibbo-suppressed D — >■ PP decays in units of 10 



Modes 


Br(diagram) 


Br(pole) 


Br(exp) 


Br (this work) 


D" ^ n^K" 


0.67±0.02 


0.6 ±0.2 




0.69 


D" ^ n-K+ 


1.12±0.05 


1.6 ±0.4 


1.48±0.07 


1.67 


D" ^ K°rj 


0.28±0.02 


0.22 ±0.05 




0.29 


D" ^ K^r,' 


0.55±0.03 


0.5 ±0.1 




0.50 


D+ -> TT+K" 


1.98±0.22 


1.7 ±0.5 




2.38 


D+ -> n''K+ 


1.59±0.15 


2.2 ±0.4 


1.72±0.19 


1.97 


D+ ^ K+r) 


0.98±0.04 


1.2 ±0.2 


1.08±0.17'" 


0.66 


D+ -> K+r,' 


0.91±0.17 


1.0±0.1 


1.76±0.22'' 


1.14 


D-^ -^ K+K° 


0.38±0.04 


0.7 ±0.4 




0.63 



"Data from [33] 
''Data from |33] 



similar to those obtained in [H [8] , and 



C2(^7r) = 1.26, C2{KK) = 1.27, 



(14) 



for the Vl^-exchange and W^-annihilation amplitudes indicate minor SU(3) breaking effects from the scale running. As 
pointed out in [TU], to account for the D° -^ n+n^ and D'^ —> K+K" branching ratios, their M^-exchange amplitudes 
must be different in both the magnitudes and strong phases. In our parametrization, the distinction in the magnitudes 
is achieved by differentiating Xq from Xs a-nd the decay constant /^r from /^^ , while the distinction in the strong phases 
is achieved by the Glauber phase S'^r- Another significant improvement appears in the prediction for the D+ — > tt+t] 
branching ratio as indicated in Table In] The branching ratios obtained in [TJ [S] are small, because the emission 
and annihilation amplitudes extracted from the Cabibbo-favored processes lead to a destructive interference in this 
Cabibbo-suppressed mode. In our framework the additional Glauber phase associated with the pion changes the 
destructive interference into a constructive one. This provides another support for the SU{3) breaking effect caused 
by the Glauber phase. 

The scale-dependent Wilson coefficients do improve the overall agreement between the theoretical predictions and 
the data involving the rj' meson. Taking the D+ — ^ n+rj' decay as an example, the corresponding parameters are 
given by ai{Trrj') = 1.12 and a-ziTrrj') = O.SOe'^"*^-^ for the emission amplitudes with the 77' meson emitted from the 
weak vertex, and Ci^'Krl') = 1.32 for the annihilation amplitudes. Compared to Eqs. (13) and (14), it is clear that 



the lower energy release involved in the D+ -^ n+rj' decay with heavier final states increases the Wilson coefficients 



Benefited from tlie scale-dependent Wilson coefficients, x^ from all the ry' involved modes is reduced by nearly 10, 
relative to x^ in the previous fits in the pole model ^. 

III. DIRECT CP ASYMMETRIES 

In this section we predict the direct CP asymmetries in the D -^ PP decays by combining the short-distance dynam- 



ics associated with the penguin operators and the long-distance parameters in Eq. (12). The direct CP asymmetry, 
defined by 

' T{D ^ f) +T{D ^ f) ^ ' 

exists only in singly Cabibbo-suppressed modes. 

A. Tree-level CP violation 

A Z? — >■ PP mode, receiving both contributions proportional to A^ and As, has the decay amplitude 

A = XdAd + \sAs. (16) 

The interference between the two terms in A leads to the tree-level direct CP asymmetry 

^ , Im(A^As) Im(^*A) ..„. 



where the relation A^-I-As « has been inserted into the denominator. Except for the final states with the j]^'^ meson, 
that contains both the qq and ss components, four other channels D^ — >■ K^K , D"*" — >■ K~^K , D+ — >■ 7r+ir°, and 
_D+ — > tt'^K'^ also exhibit the tree-level CP violation, though their values have been found to be small [10]. We can 



predict the tree-level CP violation in the D^ — )■ K^K decay as indicated in the Table 



IV 



which is not attainable in 



the diagrammatic approach in the SU{3) limit ^U\. This is one of the advantages of our formalism that has taken 
into account most of SU{3) breaking effects. The CP asymmetries in this mode, including the CP violation in the 

K'^-K mixing, would be possibly measured. 

B. Penguin- induced CP violation 

We extend our formalism to the penguin contributions to the D — > PP decays. For this purpose, the effective 
Hamiltonian for singly Cabibbo-suppressed charm decays 



^eff-^ 



E K;K,(Ci(a*)OKm) + C2(m)02'(m)) - KlKfc lj2 C.(m)0.(m) + Csgifi)OsM 

q—d,s \i—3 / 



(18) 



is considered, with the QCD penguin operators 



q'—u,d^s q'—u.d,s 

05= E {UaCa)v-A{q'l3q'fi)v+A, Oe = ^ (MaC/3)y-A (g^Qa) V-HA, (19) 

q'—u,d,s q'—u,d.s 



and the chromomagnetic penguin operator, 



87r2 



Osa = ^m.ua^^il + -f5)T-G^^''c, (20) 



T" being a color matrix. The explicit expressions for the Wilson coefficients C^^Q^fj.) in terms of the running coupling 
constant as(/i) are given in Appendix [B| 
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FIG. 2: Four topological penguin diagrams contributing to _D — >■ PP decays: (a) the penguin color-favored emission amplitude 
Pt, (b) the penguin color-suppressed emission amplitude Pc, (c) the penguin exchange amplitude Pe, and (d) the penguin 
annihilation amplitude Pa- Note that Pt{Pc) is not truly color- favored (color-suppressed), but possesses the same topology 
as T{C). 



Below we sort out the short-distance dynamics associated with the above penguin operators. The contributions 
from the {V — A){V — A) current in the operators O3.4 can be directly related to those from the tree operators by 
replacing the Wilson coefficients. The relations between the hadronic matrix elements from the {V — A){V + A) or 
{S + P){S — P) current and from the {V ~ A){V — A) current are subtler. For the penguin color-favored emission 
amplitudes, we have 



Pt = a^{^i){P2\{qq)v-An{Pi\(uc)v-A\D)+a5{lx){P2\{qq)v+A\Q){Pi\{uc)v-A\D), 
with the Wilson coefficients 

C4(m) 



(21) 



a5(M) = C5{fj.) + 






(22) 



where the minus sign in front of the Wilson coefficient a^ arises from the equality (P2I ('?'?) y+A|0) = — {P2\{qq)v -a\Q) ■ 
We write the penguin color-suppressed emission amplitude as 



Pc = a4{p){P2\iuq)v-A\0){Pi\{qc)v-A\D)~2ae{p){P2\{uq)s+p\0){Pi\ms-p\D) 

with the Wilson coefficients 

" 1 



(23) 



04 (/i) = Ci{^) + C-i{^i) 



+ Xnje 



i(f) 



Xnfe 



i(p 



The chiral factor 



2m% 



mc{mu + mg)' 



(24) 



(25) 
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with the M-quark (g-quark) mass m^ {rnq), exhibits formal suppression by a power of l/rric-, but takes a value of 
order unity actually. Note that the (5* — P){S + P) nonfactorizable amplitude cannot be simply related to the 
iy — A){V — A) nonfactorizable amplitude. However, the PQCD factorization formulas for heavy-meson decays 
reveal strong similarity between them in the dominant small parton momentum region, which supports the 



parametrization in Eq. (24 1. The chiral factor for the ri^'> meson receives an additional contribution from the axial 
anomaly, as illustrated in Appendix [C| 

The treatment of the penguin exchange amplitude Pe and the penguin annihilation amplitude P4 demands a more 
careful formulation. It is obvious that the matrix elements of the operator O3 (O4) are similar to those of O2 (Oi). 
As discussed in Sec. II, the factorizable contributions to the above matrix elements vanish in the SU{?>) limit because 
of the helicity suppression, and have been neglected. Considering only the nonfactorizable contributions, O4 gives 
rise to part of Pa- The operator Og also contributes to Pa, to which the helicity suppression applies. Therefore, its 
matrix elements can be combined with those of O4 under the equality 

{Pi[qq')P2{q'q)\{uc)v-A{q'q')v+A\D{cu)) = (Pi(<zg')P2('?'g)IMy_^(g'(7V-A|^(cS)), (26) 



which has been verified by the PQCD factorization formulas [34]. Equation (26) is attributed to the fact that only 
the vector piece V in the {q'q')v±A currents contributes to the production of two pseudoscalar mesons. Thus the 
QCD-penguin annihilation amplitude is parameterized as 

Pa = [C,[^i) + C,{ii)]xte<^fDml (^^^) • (27) 

For the amplitude Pe, the helicity suppression does not apply to the matrix elements of O^^, so the factorizable 
contributions exist. The nonfactorizable contribution from O5 vanishes in the SU{'ii) limit, as verified by the PQCD 
factorization formulas [34, and will be neglected. We then apply the Fierz transformation and the factorization 
hypothesis to 

(n(g'9l^2(uq-')|(uc)y_^(gV)y+A|^(cg))=-2(Pi(g'g)P2(V)l(^g')s+p|0)(0|(9'c)s-p|i?(cg)). (28) 

The scalar form factor involved in the matrix element {Pi{q'q)P2{uq')\{uq')s+p\Q) has been assumed to be dominated 
by scalar resonances, and fixed in 35 . It has been demonstrated that the similar formalism works for the analysis of 
the r -^ KttVt decay [3S] . Combining the nonfactorizable contribution from O3 , Pe is written as 



PE^C,{^^)xl,e<^^fornl[l^^ 



Ca(/.) ■ ""'^^^ 



Nr_ 



9sBs{ml)msfsfD^, (29) 



where gs is an effective strong coupling constant, fs and ms are the decay constant and the mass of the scalar resonance 
particle S, respectively, and the function Bg represents the Breit-Wigner propagator of the scalar resonance, 

Bs{q^) = -2 T-^. Y^Y^ . (30) 

q^ - mi. + imsTsyq^) 

In the above expression q'^ denotes the invariant mass squared of the P1P2 final state, and the strong phase from the 
width r5(g^) provides a major source to direct CP asymmetries. More detail for the treatment of the scalar matrix 
element {Pi{q' q)P2{uq')\uq'\Q) can be found in Appendix [d| 

C. Quark Loops and Magnetic Penguin 

To complete the penguin contributions to the D — > PP decays, we include the quark loops from the tree operators 
and the magnetic penguin as displayed Fig. [3j The former up to next-to-leading order in the coupling constant can 
be absorbed into the Wilson coefficients [TT] 

C,,,{n) ^ C3.5(A^)-f^ E ^C(^nM,a^)), 

q—a,s 

c,A^^) -^ c,A^^) + ^ ^ ^c^^\f,, {i')l (3i) 

QTT ^ Ab 

q—a,s 



10 




0-2 / 9 





' ■3000000^ > 
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FIG. 3: Quark-loop and magnetic-penguin contributions to D — >■ PP decays. 



with (P) being the averaged invariant mass squared of the virtual gluon emitted from the quark loop. The function 
(7(9) is written as 



c<^^\^^Al')) 



G^'^\^,,{f)) 



^2(/i), 



with the function 



, , ^ /■! ml -.t(1 - x){P) 

G(?)(^,(/2)) = -4/ dxx{l-x)\n^ ^ '—^. 

Jo f^ 



(32) 



(33) 



Because of the smallness of the associated CKM matrix elements, we have ignored the quark loops from the penguin 
operators. 



Using the unitarity relation of the CKM matrix, we reexpress the second term in Eq. (31 1 as 



"^(m) 



Stt 



A. 



Ad 
Ab 



(c'^'^^ — C^''-'^ — c'^-* 



(34) 



The term proportional to A^ contributes to the tree amplitudes, which is, however, numerically negligible, 
(as/8)7r(C*^'^^ — C'-'*-') ~ O(10~^) <^ Ci^2- It is the reason why we did not include the quark loops in the analy- 
sis of the branching ratios. Keeping only the second term in Eq. (34), Eq. (31) reduces to 



C3,5(m) ^ C3^5(M) + ^C(^nM,(^')) 



SirNr 



C4,6(m) -> C4.6(M)-^C'(^nM,(^'))- 



(35) 



We have confirmed that no much numerical difference arises from replacing C^*' by C^'^' . 

The magnetic-penguin contribution can be included into the Wilson coefficients for the penguin operators following 
the substitutions [TT] 



C3,5(M)^C3,5(M) + ^|fc-(M), 
C4,6(m) ^ ^4,6(M)-^|^Cst(M), 



(36) 



with the effective Wilson coefficient C^^ = Csg + C5. To predict the direct CP asymmetries, we adopt (P) w 
(Pi/2 + Pij'£f' ~ ?7i|)/4. This choice of (Z^) corresponds to the kinematic configuration, where the spectator in the 
decay takes half of the P\ meson momentum, consistent with the scale parametrization in Eq. (81), and a valence quark 
in the P2 meson also carries half of its momentum. We have varied (/^) between rnjj/2 and mjj/lO, and observed 
only slight change in our predictions for the direct CP asymmetries. The contribution from the magnetic penguin is 
much smaller than that from the quark loops, Cg^ ^ C^'^^{fj,, (P))- It turns out that the integral in Eq. (33) generates 
a large imaginary part as rn^ ^ (P), and a sizable shift of the strong phases in the penguin amplitudes. 
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It is seen that all the important penguin amplitudes have been formulated without introducing additional free 



parameters. Compared with [9], the similarity is that the effective Wilson coefficients in Eqs. (35 1 and (36 1 have 
been employed. The difference appears in the derivation of the involved hadronic matrix elements. Both the leading- 
power and subleading-power penguin contributions, including their strengths and strong phases, are determined in our 
formalism. The leading-power penguin contributions were calculated in QCDF, but the subleading-power ones were 
estimated in the large Nc assumption in [9^ . This is the reason why only the strengths of the penguin contributions 
were obtained in [5]. 

D. Numerical Results 

We adopt m^ = m^; — rfiq — 5.3 MeV, and Wg = 136 MeV for the quark masses at the scale 1 GeV in the numerical 
analysis. The inputs for the involved decay constants and transition form factors are taken to be the same as in 
[8j. In most modes Pc which is enhanced by the chiral factor, and Pe from the {S — P){S + P) current are the 
major QCD-penguin amplitudes. Our predictions for the direct CP asymmetries in the D — > PP decays are listed in 
Table |IV[ which differ from those derived in other approaches, and can be confronted with future data. The result 



for D'^ -^ K'^K is dominated by the tree-level CP violation, due to smallness of the involved penguin annihilation 
contribution. If considering the isospin symmetry breaking from unequal u and d quark masses, the chiral factors for 
the uu and dd components of the i:^ meson will be different. The corresponding penguin contributions then do not 
cancel exactly, and a nonzero direct CP asymmetry would appear in the D^ — > tt+tt'' decay. If turning off the Glauber 
phase in the D^ — >• tt+tt" decay, the direct CP asymmetry v4^p(7r"'"7r~) decreases from 0.58 x 10""^ to 0.49 x 10"'^. 
The 30% change in the branching ratio (see Table [TT|) and the 15% change in the direct CP asymmetry indicate that 
the Glauber phase is not a negligible mechanism. 

We predict the difference between the direct CP asymmetries in the D^ — > K^K^ and D^ -^ tt+tt" decays, 

AAcp = -1.00 X 10"^ (37) 

which is in the same sign as the data, but an order of magnitude smaller than the central value of the LHCb data in 
Eq. (II]), and lower than the CDF data in Eq. (pi). To check whether our prediction makes sense, we write A^cp as 

/\'pKK\ l-pTTTrl \ 

AAcP = -2rsin7(^[^sin<5^^ + L_isin5'^-j, (38) 

with r ~ |Ah|/|Ad| — |A6|/|As|, A, — V*gVuq, 7 being the CP violating weak phase, and S^^ {S"^^) being the relative 
strong phase between the tree amplitude f^^ fj-im-^ ^^^ ^^^ penguin amplitude V^^ (p'^'r) for the I?° — )■ K~^K~ 
(D" -^ TT+TT") mode. It is easy to find r w 0.7 x 10"^ and 7 = (73l25)° ^om the Particle Data Group 2010 [36], i.e., 
2rsin7 w 1.3 x 10^^. For |7^|/|T| ~ 0(1) and sin (5 ~ 1, AAqp could reach only a few times of 10"'^. The topological 
amplitudes for the £>" -^ tt+tt" and K^K~ modes obtained in this work are given, in unit of 10~^ GeV, by 

(39) 



(40) 
which lead to the total tree and penguin amplitudes and their ratios 

r'''' = 2.14e-*i'*°, 7''''' = 1.40e*i2i° ' = o.66e*"'*°, 

r^^'^^ = 3.94e-^i«°, T'^'^ = 1.79e*i"»^ L_ = 0.45e'"i°. (41) 

It is a concern for our formalism whether all the important nonperturbative sources in D meson decays have 
been identified. If not, the simple replacement of the Wilson coefficients may not work for estimating the penguin 





rriTTTT 


= 2.73, 


E'^'' = 0.82e" 








rpKK 


= 3.65, 


E^^ = 1.2e- 


5 




pTTTT _ 

^c - 


= 0.87e'^^^°, 


^E — 


0.81e'"i°, 


pTTTT 


= 0.25e-*43 


-,KK 

C 


= 1.21e^"5°^ 


pKK _ 
^E - 


0.87e'i"°, 


pKK 


= 0.45e-*^° 



contributions. As explained in the previous section, the outcomes in Eq. ( 12 ) are all reasonable, so most of the 
important nonperturbative sources should have been taken into account. The uncertainty of our formalism arises 
mainly from the parametrization of the scalar form factors in Pe- The difference of the direct CP asymmetries 
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TABLE IV: Direct CP asymmetries in D — >■ PP decays in the unit of (10"'^). Aq^ denotes the tree-level CP asymmetry and 
A^QP denotes the CP asymmetry arising from the interference between the total tree and penguin amplitudes. Results from the 
analysis including FSI effects gGj, and from the topological-diagram approach [TOl are presented for comparison. 



Modes 


Acp(FSl) 


j4cp (diagram) 


/I tree 


/I tot 

^CP 


D° ^ TT+TT- 


0.02±0.01 


0.86 





0.58 


D° -^ K+K- 


0.13±0.8 


-0.48 





-0.42 


D'' ^ ^V° 


-0.54±0.31 


0.85 





0.05 


D° -> if «Z" 


-0.28 ±0.16 





1.11 


1.38 


D'' ^ ^% 


1.43±0.83 


-0.16 


-0.33 


-0.29 


D« ^ ^''77' 


-0.98±0.47 


-0.01 


0.53 


1.53 


D" ^7777 


0.50±0.29 


-0.71 


0.29 


0.18 


D« ^ 7777' 


0.28 ±0.16 


0.25 


-0.30 


-0.94 


D+ ^ tt+ttO 













D+ ^ A'+l?" 


-Q.51±0.30 


-0.38 


-0.13 


-0.93 


D+ -s- 7r+77 




-0.65 


-0.54 


-0.26 


D+ ^ 7r+77' 




0.41 


0.38 


1.18 


D+ ^ ^07^+ 




0.88 


0.32 


0.39 


i3+ ^ ^+7f° 




0.52 


0.13 


0.84 


D+ ^ ^"+77 




-0.19 


0.80 


0.70 


D+ -> if +77' 




-0.41 


-0.45 


-1.60 



changes from —1.0 x 10~^ to —1.6 x 10~^, when the effective strong coupling gs = gK-'K-w is increased by a factor 
of 2. If varying the strong phases in the penguin ampHtudes arbitrarily, A^cp would be at most —1.4 x lO^'^. 
Varying the input parameters, such as the quark masses frig — 4.0 ^ 6.5 MeV and nig = 108 ^ 176 MeV, the mass 
M^^(i370) = 1200 ~ 1500 MeV and the width r/^(i37o) = 200 - 500 MeV of the dominant scalar resonance /o(1370), 

and the CKM matrix elements I Kb I = (3.89±0.44)xl0"3 ^nd 7 = ( 73^25 )°, we have A Acp = (-0.57 1.87)xl0-3. 

In conclusion, if the central values of the LHCb and CDF data persist, they can be regarded as a signal of new physics. 
At this moment, there is plenty of room for accommodating the LHCb and CDF data by means of models extended 
beyond the SM. Those models with additional sources of CP violation are all suitable candidates, like supersymmetry 
models, left-right models, models with the leptoquark, the fourth generation, or the diquark, and so on. New physics 
effects on the direct CP asymmetries in D — )■ PP decays can be explored by combining corresponding short-distance 
Wilson coefficients with the long-distance parameters determined in this work. 

IV. CONCLUSION 

In this paper we have proposed a theoretical framework for analyzing two-body nonleptonic D meson decays, based 
on the factorization of short-distance (long-distance) dynamics into Wilson coefficients (hadronic matrix elements of 
four-fermion operators). Because of the small charm quark mass just above 1 GeV, a perturbative theory for the 
hadronic matrix elements may not be valid. Our idea is to identify as complete as possible the important sources 
of nonperturbative dynamics in the decay amplitudes, and parameterize them in the framework of the factorization 
hypothesis: we have considered the evolution of the Wilson coefficients with the energy release in individual decay 
modes, which depends on the scale A of the soft degrees of freedom in the D meson. The hadronic matrix elements of 
the four-fermion operators have been parameterized into the strengths x's and the strong phases (p's. It is crucial to 
introduce the Glauber strong phase St^ associated with a pion in the nonfactorizable annihilation amplitudes, which is 
due to the unique role of the pion as a Nambu-Goldstone boson and a quark-anti-quark bound state simultaneously. 
The flavor 51/(3) symmetry breaking effects have been taken into account in the above framework appropriately. 
Fitting our parametrization to the abundant data of the D — >■ PP branching ratios, all the nonperturbative parameters 
have been determined. 

It has been shown that our framework greatly improves the global fit to the measured D — >■ PP branching ratios: 
the evolution of the Wilson coefficients improves the overall agreement with the data involving the 77' meson, and 
the Glauber phase resolves the long-standing puzzle from the D'^ — > tt+tt^ and D'^ — > K^K^ branching ratios. The 
Glauber phase has also enhanced the predicted D^ — >■ 7r+77 branching ratio, giving much better consistency with the 
data. The value of A is in the correct order of magnitude for characterizing the soft degrees of freedom in the D meson. 
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and the Glauber phase S'jr agrees with that extracted from the data for the direct CP asymmetries in the B — >■ ttK 
decays. The similarity of the x values for the VF-exchange and for the W-annihilation conforms our expectation, 
because these two topologies can be converted to each other via the Fierz transformation. The difference between the 
X values for the strongly produced light-quark and strange-quark pairs implies the significant SU{3) breaking effects 
in the D — ^ PP decays. 

Once having determined the nonperturbative parameters from the fit to the branching ratios, the replacement of 
the Wilson coefficients works for estimating the penguin contributions. For those penguin amplitudes, which cannot 
be related to tree amplitudes through the above replacement, we have shown that they are either factorizable or 
suppressed by the helicity conservation. If they are factorizable, such as the scalar penguin annihilation contribution, 
data from other processes can be used for their determination. We are then able to predict the direct CP asymmetries 
in D — >■ PP decays without ambiguity. It has been found that the strong phases from the quark loops and the scalar 
penguin annihilation dominate the direct CP asymmetries. Many of our predictions in Table |IV| can be confronted 
with future data. In particular, we have predicted AAcp = —1.00 x 10~^, which is lower than the LHCb and CDF 
data. As pointed out at the end of the previous section, the uncertainty of our formalism arises mainly from the 
parametrization of the scalar form factors in Pe- However, even increasing the associated effective strong coupling 
by a factor or 2, AAcp remains of order 10"'^ . We conclude that the LHCb and CDF data will reveal a new-physics 
signal, if their central values persist. 
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Appendix A: Amplitudes 

Different decay modes receive contributions from different topological amplitudes. Take the D" — > tt+tt^ and 
D'^ -^ K^K^ modes as examples, whose expressions are give by 



A{D' -> vr+.-) = ^ 



A{D^ ^K+K-) = 



Xd (T + E)- X, (Pc + 2P'i + Pi 






JO 



Gp 
Gp 



Xs (T + E)- Xh (Pc + P1 + P1 + Pk) 



v:,Vus 



a,i^,){ml - m],)fKF^''(m],) + C2(M)e<xf /i^m^ 



f2 
2 lK_ 

q jiy-D n 



ap^{^x){ml-ml)jKF^^{m\)+ap,{p)e<XqlDnil^-^ 



+apAf^)e'^^xffDml^+C^{f,)e^^fxffDml^ 






(Al) 



with the coefBcients 
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aPcip) = [a4(M)+a6(M)rJ, 

a',{^,) = Cdp) + ^. (A2) 



Notice the Glauber phase factors e^ °^ associated with the nonfactorizable amphtudes Xq and Xq in Eq. (All 



Appendix B: Wilson coefficients 

In this Appendix we present the evolution of the Wilson coefficients in the scale jjl < nic [37] , 

Ci(Ai) = 0.2334ai-444 + 0M59a°-^™ - lM3a°"^^^ + OMila-"-^^^^, (Bl) 

C2i^i) = -0.2334ai-444 + 0.0459a°-^"^« + 1.313a°-4444 + 0.3041a-°-2222, (B2) 

C3i^i) = 0.0496a-"-2977 - 0.0608a-°-2222 + 0.0025a-°"96 _ i_23a°i''^ 
+1.313a°-4444 - 0.0184a°™23 _ 0.0076a°-^^^« - 0.0188a°-»°25 

+0.0018a"-««°4 + 0.2362ai-4i^ - 0.1986ai-444 + 0.0004ai-8"2, (B3) 

C4(m) = 0.075a-°-2977 - 0.0608«-"-2222 ^ o.0012a"0"96 + l.i79a"-4i^3 
-1.313a"-4444 + 0.0175a°™23 _ 0.014a°-^^^8 _^ o.0267a°-8025 

+0.0006a°««°4 - 0.1807ai-4i^ + 0.129a^-^^^ - 0m2a^'^"^, (B4) 

C^ifi) = -0.0252a-0-2977 ^ o.0227a-""96 + 0.0309a°-4i^3 ^ o.0251a"™23 
+0.0016a°-^^^« - 0.0045a0«°25 „ 0.0058a"«804 _ o.0338ai-4i^ 

+0.0348ai-*4* - 0.026ai-8°2^ (B5) 

C6(Ai) = 0.019a-0-29^7 _ o.0081a-""9f' + 0.0821a°-4i73 _ o.0624aO-^023 
-0.0048a°-^^^« - 0.1241a0-8025 „ o.0038a°-8804 + 0.0979ai-4i^ 

-0.1045ai-*^'' - 0.0486al■*°^ (B6) 

in terms of the running coupling constant 



with the coefficients 



47r 


1 /3ilnln(MVAy 

[ PI i^(MVAy 


/3oln(MVAy 


33-2/ 

- 3 , /3i 


-102-fiV/. 



(B7) 



(B8) 



For the active flavor number Nf — 3 and the QCD scale Aj^^ = A^ — 375 MeV, we have as{mc) — 0.407. The 
initial conditions of the Wilson coefficients at /i = nic are given, in the NDR scheme, by 

Clime) = -0.43, C2(m,)-1.22, 

Csiruc) = 0.018, Ciirric) = -0.046, (B9) 

C5(mc) = 0.013, Ceim,) = -0.044. 

Another useful set of values contains 

Ci(lGeV) = -0.51, C2(lGeV) = 1.27, 

C3(l GeV) = 0.028, ^(l GeV) = -0.065, (BIO) 

C5(lGeV) = 0.021, 6*6(1 GeV) = -0.058, 

which are evaluated at the typical energy release in two-body nonleptonic D meson decays. 

The effective Wilson coefficient C|g for the magnetic penguin operator is known only to the leading-logarithm 
accuracy. Because of its smallness, we shall not consider the evolution of Cgg (/i), but approximate it by the constant 

Ctim,) = -0.11. (Bll) 
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Appendix C: Chiral factor 

The chiral factor in Eq. (l25) comes from the relation between the matrix elements of the scalar (pseudoscalar) 
current and of the vector (axial vector) current. Using the equations of motion, we derive 



-idt'{P2\ui^l5q\0) _ . fP2m 



2 
2 '"Pa 



{P,\u^,q\0) = ^^^^2"'^^" = -^^^V^, (CI) 

rriu + rUq rriu + rriq 

{Pi\qc\D) ^ = ^0 (™P2). (C2) 

mc rric 

where the light quark mass has been neglected in the second expression. 

For the rj^'' meson, the chiral factor receives an additional contribution from the axial anomaly in the following 
equation of motion |38j 

5^(g7,.75g) = 2m,{qzj,q) + ^GG, (C3) 

in which G represents the giuon field tensor, G is its dual, and q denotes the quark fields u, d, and s. In the present 
analysis only the u quark contributes to the amplitude Pc- We express the physical states rj and rj' as the linear 
combinations of the flavor states rjq and rjg via a mixing angle </>, 

^ / cos (/> - sin (/) \ ( Vq\ /(-;4^ 

\ sin (j) cos (p I \ rjs I 

The KLOE collaboration has extracted the value (j> — (40.4 ± 0.6)° recently [32], which is consistent with the phe- 
nomenological result in |38l. 




Equations ( C3 ) and ( C4 1 then lead to the matrix elements 



(77|m75«|0) = ^(^-^m2/,cos</>-cos</>(r;,|^GG|O)+sin0(,7,|^GG|O)^, (C5) 

W\mj5u\0) - ^('-^m2,/<,sin0-sin0(r;,|^GG|O)-cos0(7,,|^GG|O)V (C6) 



with the rjq decay constant fq. For the expressions of the anomalies 



{Vq\"^GG\0) ^ V2fqa', {rj,\"^GG\0)=yfqa', (C7) 



I^GGIO) = V2fqa\ (,?,|£ 
the involved parameters have been determined to be fq — 1.07 f-^, a^ — 0.265 GeV^, and y = 0.81 



Appendix D: Scalar matrix elements 

In this Appendix we fix the scalar matrix element {PiP2\qiQ2\0) appearing in Sec. Ill C, which is dominated by the 
contribution from the lowest scalar resonance |35j . 

{PlP2\qiq2\0) = {PiP2\S){S\qiq2\0) = gsBsiq^)msfs, (Dl) 



as shown in Eq. (30). The scalar decay constant fs is defined via (S'|g29i|0) = msfs [lO]- Many scalar mesons have 
been observed. It is still controversial whether the lightest scalar nonet with the mass smaller than or close to 1 GeV 
are primarily the four-quark or two-quark bound states. Since k, ao(980), and /o(980) may be the four-quark states, 
we consider the scalar resonances ao(1450) for the ■n-q, -kt]' , and [KKY^ final states, /o(1370), /o(1500), and /o(1710) 
for the 7r°7r°, tt+tt", K^K^, K+K', and rjrj'-'^ final states, and ii:o(1430) for the Ktt and K-q'^''^ final states due to the 
symmetries of strong interaction. /o(1370), /o(1500), and /o(1710) are the mixtures of nh = {uu + d(J)/-\/2, ss and 
glueball, but there still exists controversy on the components of each scalar meson. We employ the mixing matrix 
obtained in |41j in this work. Note that there are no s-wave isospin-1 resonances for the tt tt'^ system under the 
isospin symmetry. 

We derive (i^7r|su|0) as an example. Inserting ii'Q(1430) as the intermediate resonance, the matrix element is 
written as 

{Ktt\su\0) = (i^TTlif* (1430)) (i^o* (1430) |5-u|0) = gK*K.BK:M^)mK^jK'^, (D2) 
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where the effective strong couphng constant qk'Ktt — 3.8 GeV can be obtained directly from the measured Kq (1430) ^• 
Ktt decay 02], and Bk' is the iiTg (1430) propagator in the Breit-Wigner form [35] 

Bk'M^) = -^ ^— ^ T^^^' (D3) 

with the momentum-dependent width 

„2n _ T^__ n__ ^„2 






2, _ mK^ PK{q^) 

PK{q') = -^^\'l\q\ml.m\\ 
A(a;, y, z) = x^ + y^ + z^ - 2xy - 2yz - 2xz. (D4) 

The values of the mass rriK* and the decay width Tx' for the scalar resonance iirQ(1430) are taken from the Particle 
Data Group 2010 131]. 

The relevant scalar decay constants are given, in QCD sum rules, by j40) 

/„fi(lGeV) = (460±50)MeV, 
/s5(lGeV) = (490±50)MeV, 
/k* (1430) (1 GeV) = (445±50)MeV, (D5) 

where fnn is the scalar decay constant for ao(1450) and for the nn components of the three iso-singlet scalars, and 
fss for the ss components of the three iso-singlet scalars. We assume that the quark, but not glueball, components of 



scalars dominate the D meson decays. With lack of experimental data, the effective coupling constant gs in Eq. (Dl ) 
is approximated by gs — gK'KTr in the flavor SU{S) limit. 
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